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Abstract 

Combining work of Serre and Zumbrun, Benzoni-Gavage, Serre, and Zumbrun, and 
Texier and Zumbrun, we propose as a mechanism for the onset of cellular instability 
of viscous shock and detonation waves in a finite-cross-section duct the violation of 
the refined planar stability condition of Zumbrun-Serre, a viscous correction of the 
inviscid planar stability condition of Majda. More precisely, we show for a model 
problem involving flow in a rectangular duct with artificial periodic boundary conditions 
that transition to multidimensional instability through violation of the refined stability 
condition of planar viscous shock waves on the whole space generically implies for a 
duct of sufficiently large cross-section a cascade of Hopf bifurcations involving more and 
more complicated cellular instabilities. The refined condition is numerically calculable 
as described in Benzoni-Gavage-Serre-Zumbrun. 



1 Introduction 

It is well known both experimentally and numerically [BEl iMTl IBMRl iFWj iMTl lAYT\ 
\AT\ \F1\ \F2\ IKSj that shock and detonation waves propagating in a finite cross-section 
duct can exhibit time-oscillatory or "cellular" instabilities, in which the initially nearly 
planar shock takes on nontrivial transverse geometry. Majda et al [MRU IMR21 lAMj have 
studied the onset of such instabilities by weakly nonlinear optics expansion of the associated 
planar inviscid shock in the whole space. More recently, Kasimov-Stewart |KS| and Texier- 
Zumbrun [TZ21 ITZ3[ ITZ3| have studied these instabilities as Hopf bifurcations of flow in 
a finite-cross-section duct, associated with passage across the imaginary axis of eigenvalues 
of the linearized operator about the wave. 

In this paper, combining the analyses of [BSZl[TZ4llZ4] . we make an explicit connection 
between stability of planar shocks on the whole space, and Hopf bifurcation in a finite cross- 
section duct, by a mechanism different from that investigated by Majda et al. Specifically, 
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we point out that violation of the refined stability condition of \ZS\ \Z1\ IBSZ| . a viscous 
correction of the inviscid planar stability condition of Majda |Mlj - [M4] . is generically as- 
sociated with Hopf bifurcation in a finite cross-section duct corresponding to the observed 
cellular instability, for cross-section M sufficiently large. Indeed, we show more, that this 
is associated with a cascade of bifurcations to higher and higher wave numbers and more 
and more complicated solutions, with features on finer and finer length/time scales. 

1.1 Equations and assumptions 

Consider a planar viscous shock solution 

(1.1) u{x,t) = u{xi — st) 
of a two-dimensional system of viscous conservation laws 

(1.2) ^t + Y^ fKu)xj = A^ii, u G M", X G M^ t € M+ 

on the whole space. This may be viewed alternatively as a planar traveling-wave solution 
of ()1.2p on an infinite channel 

C:={x: (xi,X2) G x [-M,M]} 

under periodic boundary conditions 

(1.3) n(xi,M) = n(xi,-M). 

We take this as a simplified mathematical model for compressible flow in a duct, in 
which we have neglected boundary-layer phenomena along the wall dVt in order to isolate 
the oscillatory phenomena of our main interest. 

Following [TZ2j . consider a one-parameter family of standing planar viscous shock so- 
lutions u^{xi) of a smoothly- varying family of conservation laws 

2 

(1.4) ut = He,u):=A^u-Y,F^{e,u)^^, n G 

i=i 

in a fixed channel C, with periodic boundary conditions (typically, shifts {e,u)x '■= 
Yl,f^^'^)x — s{e)Ux^ of a single equation (jl.2p written in coordinates xi — > xi — s{e)t 
moving with traveling- wave solutions of varying speeds s(e)), with linearized operators 
u=u^ ■ Profiles satisfy the standing- wave ODE 

(1.5) u' = F^{e,u) - F\e,u.). 
Let 

(1.6) Alie):= liin F^ie,^^. 

Following [ZH ITZ21 IZ4j . we make the assumptions: 
(HO) F^ €C\k> 2. 
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(HI) a{A^{e)) real, distinct, and nonzero, and ^j^±i^)) '^sal and semisimple 
for ^ G M'^'. 

For most of our results, we require also: 

(H2) Considered as connecting orbits of (jl.Sp . are transverse and unique up to 
translation, with dimensions of the stable subpace S{A\_) and the unstable subspace U{A^) 
summing for each e to n + 1. 

(H3) det(rj~, . . . , r~ ,r'^_^_i, . . . , r+, n+— u_) / 0, where rj~, . . . , r~ are eigenvectors of 
Al_ associated with negative eigenvalues and r~^^, . . . ,r~ are eigenvectors of associated 
with positive eigenvalues. 

Hypothesis (H2) asserts in particular that is of standard Lax type, meaning that the 
axial hyperbolic convection matrices A^{£) and ^L(e) at plus and minus spatial infinity 
have, respectively, n — c positive and c — 1 negative real eigenvalues for 1 < c < n, where c is 
the characteristic family associated with the shock: in other words, there are precisely n — 1 
outgoing hyperbolic characteristics in the far field. Hypothesis (H3) may be recognized 
as the Liu-Majda condition corresponding to one-dimensional stability of the associated 
inviscid shock. In the present, viscous, context, this, together with transverality, (H2), 
plays the role of a spectral nondegeneracy condition corresponding in a generalized sense 
[ZHl IZlj to simplicity of the embedded zero eigenvalue associated with eigenfunction dx^u 
and translational invariance. 

1.2 Stability conditions 

Our first set of results, generalizing the one-dimensional analysis of |Z4| . characterize stabil- 
ity/instability of waves in terms of the spectrum of the linearized operator Fixing 
e, we suppress the parameter e. We start with the routine observation that the semilinear 
parabolic equation ()1.2p has a center-stable manifold about the equilibrium solution u. 

Proposition 1.1. Under assumptions (HO)-(Hl), there exists in an neighborhood of the 
set of translates of u a codimension-p translation invariant C'^ (with respect to H^) center 
stable manifold A4cs, tangent at u to the center stable subspace Tics of L, that is (locally) 
invariant under the forward time- evolution of (11.20 ~ (ll.3p and contains all solutions that 
remain bounded and sufficiently close to a translate of u in forward time, where p is the 
(necessarily finite) number of unstable, i.e., positive real part, eigenvalues of L. 

Proof. By standard considerations [Hel ITZ2] . L{e) possesses no essential spectrum and at 
most a finite set of positive real part eigenvalues on 3?A > 0. With this observation, the 
result follows word-for-word by the argument of |Z4] in the one-dimensional case, which 
depends only on the properties of L as a sectorial second-order elliptic operator, and on 
semilinearity and translation-invariance of the underlying equations (jl.2p . □ 

Introduce now the nonbifurcation condition: 

(Dl) L has no nonzero imaginary eigenvalues. 
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As discussed above, (H2)-(H3) correspond to a generalized notion of simplicity of the 
embedded eigenvalue A = of L. Thus, (Dl) together with (H2)-(H3) correspond to the 
assumption that there are no additional (usual or generalized) eigenvalues on the imaginary 
axis other than the translational eigenvalue at A = 0; that is, the shock is not in transition 
between different degrees of stability, but has stability properties that are insensitive to 
small variations in parameters. 

Theorem 1.2. Under (H0)-(H3) and (Dl), u is nonlinearly orhitally stable as a solution 
of ()1.2p - (jl.3p under sufficiently small perturbations in fl lying on the codimension 
p center stable manifold Aics of u and its translates, where p is the number of unstable 
eigenvalues of L, in the sense that, for some a(-), all L'^, 

a{t))\LP < C{1 + ty^'^^^p'^\u{x,0) - u{x)\LinH^, 

a{t))\H2 < C{l + ty^\u{x,0) -u{x)\Linm, 

a{t) < C{l + t)-^\u{x,0) -u{x)\LinH2, 

a{t) < C\u{x,0) - u{x)\LinH2- 

Moreover, it is orbitally unstable with respect to small perturbations not lying in Mcs, 
in the sense that the corresponding solution leaves a fixed-radius neighborhood of the set of 
translates of u in finite time. 

Remark 1.3. Theorem II .21 includes in passing the result that existence of unstable eigen- 
values implies nonlinear instability, hence completey characterizes stability/instability of 
waves under the nondegeneracy condition (Dl). The rates of decay ()1.7p are exactly those 
of the one-dimensional case [Z4j . 

1.3 Bifurcation conditions 

We next recall the following result from |TZ2t ITZ Sj characterizing Hopf bifurcation of 
in terms of conditions on the spectrum of L{e). Define the Hopf bifurcation condition: 

(D2) Outside the essential spectrum of L{e), for e and 5 > sufficiently small, the 
only eigenvalues of L{e) with real part of absolute value less than 6 are a crossing conjugate 
pair l±{e) := 7(e) ± ir(e) of L(e), with 7(0) = 0, 9^7(0) > 0, and r(0) / 0. 

Proposition 1.1 ( |TZ21 ITZ3| ). Let , (|1.4p be a family of traveling-waves and systems 
satisfying assumptions (H0)-(H3) and (D2), and tj > sufficiently small. Then, for a > 
sufficiently small and C > sufficiently large, there are functions e{a), e(0) = 0, and 
T*{a), r*(0) = 27r/r(0), and a family of solutions n"(xi,t) of (fLlj) with e = e{a), 
time-periodic of period T*[a), such that 

(1.8) C^^a < sup e''l^il |n"(x,t) -u"('')(xi)| < C7a for ah t > 0. 

Up to fixed translations in x, t, for e sufficiently small, these are the only nearby solutions 
as measured in norm := 11(1 + |2^i[)/(2;)||l°°(2:) ^^'^^ '^'"c time-periodic with period 



(1.7) 



\u{x, t) — u{x 
\u{x, t) — u{x - 
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T E [To,Ti], for any fixed < Tq < Ti < +00. Indeed, they are the only nearby solutions 
of form u"'{x,t) = u"(x — a"'t,t) with u"^ periodic in its second argument. 

Proof. This result was established in Theorem 1.4, |TZ2j with (jl.Sp replaced by 

(1.9) C-^a< sup {l + \xi\)\u''{x,t) -u'^^^xi)] < Ca 

and under the further assumption that there are no eigenvalues of L{£) with strictly pos- 
itive real part other than possibly A±(e). As L(e) by standard considerations |Hel ITZ2] 
possesses at most a finite set of positive real part eigenvalues, an examination of the proof 
shows that the more general case follows by essentially the same argument, reducing by the 
Lyapunov-Schmidt reduction described in |TZ2] to a finite-dimensional equation on the di- 
rect sum of the oscillatory eigenspace associated with X± and the unstable eigenspace of L, 
then appealing to standard, finite-dimensional theory to conclude the appearance of Hopf 
bifurcation with bound (jl.Op . The stronger result of exponential localization, (jl.Sp . may be 
obtained by combining the argument of [TZ2j with the strengthened cancellation estimates 
of Proposition 2.5 [ TZ3] . As the distinction between (jl.Sp and (jl.9p is not important for 
the present discussion, we omit the (straighforward) details. □ 

Remark 1.4. Together with Theorem 11.2^ Proposition 11.11 implies that, under the Hopf 
bifurcation assumption (D2) together with the further assumption that L{£) have no strictly 
positive real part eigenvalues other than possibly X±, waves are linearly and nonlinearly 
stable for e < and unstable for e > 0, with bifurcation/exchange of stability at e = 0. 

1.4 Longitudinal vs. transverse bifurcation 

The analysis of |TZ2| in fact gives slightly more information. Denote by 

(1.10) U^f:=Y,^^j{xWjJ) 

the L(e)-invariant projection onto the oscillatory eigenspace := Span{0^}, where are 
the eigenfunctions associated with A-i-(e). Then, we have the following result, proved but 
not explicitly stated in |TZ2j . 

Proposition 1.5. Under the assumptions of Proposition [Ll\ also 

(1.11) sup e''l^il|'u" -u-n^('u<' -u)| < Ca^ for ah t > 0. 

XI 

Proof. The weaker bound 

(1.12) supil + IxiDlu" -u-U'iu" -u)\<Ca'^ for alH > 0, 

XI 

is established in the course of the Lyapunov reduction of |TZ2| : see (2.17), Proposition 2.9, 
case u; = 0. The stronger version (jl.lip follows by the same argument together with the 
strengthened cancellation estimates of Proposition 2.5 [TZ3j . □ 
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Bounds (jl.Sp and (jl.lip together yield the standard finite-dimensional property that 
bifurcating solutions lie to quadratic order in the direction of the oscillatory eigenspace of 
L{e). From this, we may draw the following additional conclusions about the structure of 
bifurcating waves. By separation of variables, and X2-independence of the coefficients of 
L{e), we have that the eigenfunctions of L{e) decompose into families 

(1-13) e'^^-tPix,), e = ^, 

associated with different integers k, where M is cross-sectional width. Thus, there are two 
very different cases: (i) (longitudinal instability) the bifurcating eigenvalues A±(e) are asso- 
ciated with wave-number /c = 0, or (ii) (transverse instability) the bifurcating eigenvalues 
A±(e:) are associated with wave-numbers ±k ^ 0. 

Corollary 1.6. Under the assumptions of Proposition u"" depend nontrivially on X2 
if and only if the bifurcating eigenvalues X± are associated with transverse wave-numbers 
±k / 0. 

Proof. For k ^ 0, the result follows by the fact that, by (jl.Sp and (jl.lip . Il{u'^ — u) is the 
dominant part of — n, and the fact that 11/ by inspection depends nontrivially on X2 
whenever 11/ ^ 0. For A; = 0, the result follows by uniqueness, and the fact that, restricted 
to the one-dimensional case, the same argument yields a bifurcating solution depending 
only on xi. □ 

Bifurcation through longitudinal instability corresponds to "galloping" or "pulsating" in- 
stabilities described in detonation literature, while symmetry-breaking bifurcation through 
transverse instability corresponds to "cellular" instabilities introducing nontrivial transverse 
geometry to the structure of the propagating wave. 

1.5 The refined stability condition and bifurcation 

Longitudinal or "galloping" bifurcation, though almost certainly occurring for detonations 
(see |TZ2[[TZ4] and references therein), has up to now not been observed for shock waves as 
far as we know (though we see no reason why they should not in general be possible) , nor has 
there been proposed any specific mechanism by which this might occur. The main purpose 
of the present paper, as we now describe, is to point out that for transverse or "cellular" 
bifurcations, to the contrary, there is a simple and natural mathematical mechanism, closely 
related to the inviscid stability theory for shocks in the whole space, by which they can and 
likely do occur. 

1.5.1 The inviscid stability condition 

Inviscid stability analysis for shocks in the whole space centers about the Lopatinski deter- 
minant 

(1.14) A(e,A):=(7^^ ••• 7e++i ••• n+ AM + iel/']), 
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M , X = 'J + ir C, T > 0, a spectral determinant whose zeroes correspond to 
normal modes e^''e^^^'^w{xi) of the constant-coefficient linearized equations about the dis- 
continuous shock solution. Here, {TZpj^i, . . . , Tin} ^'^^d {Tl-i, • • • , ^p_i} denote bases for the 
unstable/resp. stable subspaces of 

(1.15) A+ii, A) := (A/ + i^dfiu±))idf\u±))-\ 

Weak stability |A[ > for r > is clearly necessary for linearized stability, while strong, 
or uniform stability, |A|/|(^,A)| > cq > 0, is sufficient for nonlinear stability. Between 
strong instability, or failure of weak stability, and strong stability, there lies a region of 
neutral stability corresponding to the appearance of surface waves propagating along the 
shock front, for which A is nonvanishing for 5RA > but has one or more roots {S,o,Xo) 
with Ao = To pure imaginary. This region of neutral inviscid stability typically occupies an 
open set in physical parameter space |MH IM2t \M3\ |BRSZ', 'Zlj IZ2j . For details, see, e.g., 
[ETl IHH IMl IMHI [Mil Eiil EM EiHl EHl EH ESI El, BRSZ] , and references therein. 

It has been suggested [MRU IMR21 lAMj that nonlinear hyperbolic evolution of surface 
waves in the region of neutral linear stability might explain the onset of complex behavior 
such as Mach stem formation/kinking of the shock. We pursue here a variant of this idea 
based instead on interaction between neglected viscous effects and transverse spatial scales. 

1.5.2 The refined stability condition 

Viscous stability analysis for shocks in the whole space centers about the Evans function 

D{1 A), 

^ G M^, A = 7 -|- € C, r > 0, a spectral determinant analogous to the Lopatinski 
determinant of the inviscid theory, whose zeroes correspond to normal modes e'^^e^^^'^w{xi), 
of the linearized equations about u (now variable-coefficient), or spectra of the linearized 
operator about the wave. The main result of [ZS] , establishing a rigorous relation between 
viscous and inviscid stability, was the asymptotic expansion 

(1.16) Z)(|,A) = 7A(e,A) + o(|(|,A)|) 

of D about the origin {(,,X) = (0, 0), where 7 is a constant measuring tranversality of u as a 
connecting orbit of the traveling-wave ODE. Equivalently, considering -D(^, A) = D(p£^o, pXo) 
as a function of polar coordinates {p, ^O) Aq), we have 

(1.17) D\p=o = and {d/dp)\p=oD = ^Ai^o, Ao). 

An important consequence of ()1.16p is that weak inviscid stability, |A| > 0, is necessary 
for weak viscous stability, \D\ > (an evident necessary condition for linearized viscous 
stability). For, (I1.16|) implies that the zero set of D is tangent at the origin to the cone 
{A = 0} (recall, (|1.14p . that A is homogeneous, degree one), hence enters {r > 0} if {A = 0} 
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does. Moreover, in case of neutral inviscid stability A(^0;^'^o) = 0, (^0)'^''"o) (OiO)> one 
may extract a further, refined stability condition 

(1.18) P := -Dpp/D,x\p=o > 

necessary for weak viscous stability. For, (jl.l7p then impHes Dp\p=Q = 7A(,^Oi^'^o) = 0, 
whence Taylor expansion of D yields that the zero level set of D is concave or convex 
toward r > according as the sign of f3; see |ZS| for details. As discussed in \ZS\ IZlj . 
the constant P has a heuristic interpretation as an effective diffusion coefficient for surface 
waves moving along the front. 

As shown in |ZSl IBSZj . the formula ()1.18p is well-defined whenever A is analytic 
at (Coj^'T"); in which case D considered as a function of polar coordinates is analytic at 
(0, ^0) ^''"o)) and iro is a simple root of A(^0;')- The determinant A in turn is analytic 
at {Co, iro), for all except a finite set of branch singularities tq = ^o^j- As discussed in 
|BSZ1 IZ21 IZ3j , the apparently nongeneric behavior that the family of holomorphic functions 
A"^ associated with shocks (ti^,nl) have roots (^Q,iro(e)) with iro pure imaginary on an 
open set of e is explained by the fact that, on certain components of the complement on the 
imaginary axis of this finite set of branch singularities, A^(,^q, •) takes the imaginary axis 
to itself. Thus, zeros of odd multiplicity persist on the imaginary axis, by consideration of 
the topological degree of A^ as a map from the imaginary axis to itself. 

Moreover, the same topological considerations show that a simple imaginary root of 
this type can only enter or leave the imaginary axis at a branch singularity of A^(^'^, •) or 
at infinity, which greatly aids in the computation of transition points for inviscid stability 
|BSZl[ZTllZ2UZ3j . As described in |Z2 1 IZ3 irSelj . escape to infinity is always associated with 
transition to strong instability. Indeed, using real homogeneity of A, we may rescale by |A| 
to find in the limit as |A| oo that = |Ao|"^A(|o, Aq) = A(|o/|Ao|, Ao/|Ao|) A(0,i), 
which, by the complex homogeneity A(0, A) = AA(0, 1) of the one-dimensional Lopatinski 
determinant A(0, •), yields one- dimensional instability A(0, 1) = 0. As described in [Zl], 
Section 6.2, this is associated not with surface waves, but the more dramatic phenomenon 
of wave-splitting, in which the axial structure of the front bifurcates from a single shock to 
a more complicated multi-wave Riemann pattern. 

Example 1.2. For gas dynamnics, complex symmetry, A(^, A) = A(— ^, A), and rotational 
invariance, A(4,A) = A(— ^, A), imply that 

(1.19) A(f,ir) = A(|eT,|r|2). 

Explicit computation \Erl\ \M1\ \Z1^ yields that A(,^Oi ■) has a pair of branch points of square- 
root type, located at |toP = |.^oP(-^^ ~ 1)? where M is the downstream Mach number (? ju^ , 
where c is sound speed and u the axial particle velocity of the shock on the downstream side, 
defined as the side in the direction of particle velocity. Transition from strong stability to 
neutral stability occurs through a pair of simple imaginary zeros entering the imaginary axis 
at the branch points, and transition from neutral stability to strong instability occurs through 
escape of these zeros to infinity, with associated one- dimensional instability/wave-splitting. 
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Remark 1.7. We note in passing that one-dimensional inviscid stability A^(0, 1) 7^ is 
equivalent to (H3) through the relation 



(1.20) A^(0,A) = Adet(r]-,, 



1.5.3 Transverse bifurcation of flow in a duct 

We now make an elementary observation connecting cellular bifurcation of flow in a duct 
to stability of shocks in the whole space: specifically, to violation of the refined stability 
condition. Assume for the family the stability conditions: 

(Bl) For e sufficiently small, the inviscid shock {u'^,u'L) is weakly stable; more 
precisely, A^(1,A) has no roots > but a single simple pure imaginary root A(e) = 
ir*(e) 7^ lying away from the singularities of A^. 

(B2) The refined stability coefficient /3(e) defined in (fLTSD satisfies 3?/3(0) = 0, 
d,^p{0) < 0. 



Lemma 1.8 ( |ZSl IZT] ). Assuming (H0)-(H2) and (Bl), fore, ^ sufficiently small, there 
exist a smooth family of roots (C)'^*(C)) ^/-^(C)'^) with 

(1.21) K[0 = ^ir*{e) - + 5{e)~e + r{e,i)i\ r G C\e,0- 

Moreover, these are the unique roots of D satisfying KA > — |^|/C for some C > and 
p= |(CjA)| sufficiently small. 

Proof. This follows by the Implicit Function Theorem applied to the function D{p, Aq) := 
p^^D{p^Q, pXo), about the values {p, Aq) = (0, i^oT^{e), where ^0 is without loss of generality 
held fixed, using the facts that D expressed in polar coordinates {p, ^0, Aq) satisfies D\p=Q = 
and dpD\p=Q = A, so that Dp, Dw, and Dx all vanish at {0,S,o,S,oiT^{e)). For details, see 
the proof of Theorem 3.7, [Zl| . □ 



Corollary 1.9. Assuming (H0)-(H2) and (B1)-(B2), for e, ^ sufficiently small, there is a 
unique function £{6.) 7^ 0, £{0) = O, such that 5RAJ(^) = for e = £{C)- the generic 
case 3fJ5(0) 7^ 0, moreover, 

(1.22) £ii) ~ {m{o)/dem)l 

Proof. As a consequence of (|1.2ip . we have for some smooth G 

(1.23) ^(^) = -me)+iGie,i), 



G := {6{e) + r(e,0), whence the equation = = -^(3{e) + ^G(^,e) has a 

unique root e = £{S,) by assumption (B2) and standard scalar bifurcation theory. From 
G = 5(e)| + 0(|f|2), we find, in the generic case m{0) / 0, that d^{iG)\^^^^QQ = m{0) + 0, 
yielding (fr22]l . ' ' □ 
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Remark 1.10. As described further in [BSZ| . the above resuhs on the refined stabihty 
condition apply also in the case of "real" or partial viscosity, in particular to the physical 
Navier-Stokes equations of compressible gas dynamics and MHD. 

To (Bl) and (B2), adjoin now the additional assumptions: 
(B3) (5(0) / 0. 

(B4) At e = 0, the Evans function D{S^,X) has no roots ^ G M, 9?A > outside a 
sufficiently small ball about the origin. 

Then, we have the following main result. 
Theorem 1.11. Assuming (H0)-(H2) and (B1)-(B4), for 

^max > sufficiently small and 
each cross-sectional width M sufficiently large, there is a finite sequence < £i{M) < • • • < 
efc(M) < ••• < £max, with £k{M) ~ {^6{0)/deP{0))j^, such that, as e crosses successive 
Ek from the left, there occur a series of transverse (i.e., "cellular") Hopf bifurcations of 
associated with wave-numbers ±k, with successively smaller periods T/c(e) ''"*(0)^^. 

Proof By (B4), for |e| < 

^max sufficiently small, we have by continuity that there exist no 
roots of A) for KA > —1/C, C > 0, outside a small ball about the origin. By Lemma 
11.81 within this small ball, there are no roots other than possibly (^, A'^(.^)) with 3?A > 0: in 
particular, no nonzero purely imaginary spectra are possible other than at values X^{S,) for 
operator L{e) acting on functions on the whole space. 

Considering L instead as an operator acting on functions on the channel C := {x : 
(xi,X2) € M"*^ X [— Af, M]}, we find by discrete Fourier transform/separation of variables 
that its spectra are exactly the zeros of D{^k, A), as = ^ runs through all integer wave- 
numbers k; see ()1.13p . Applying Corollary 11.91 and using (B3), we find, therefore, that pure 
imaginary eigenvalues of L{e) with |e| < £max sufficiently small occur precisely at values 
£ = £k, and consist of crossing conjugate pairs A^(e) associated with wave-numbers ±/c, 
satisfying Hopf bifurcation condition (D2) with 

9Ai(e)~n(0))^. 

Applying Proposition II. H we obtain the result. □ 

Remark 1.12. As evidenced by decreasing periods T^, this phenomenon of increasing- 
complexity solutions is completely different from the more familiar one of period-doubling. 

Remark 1.13. Lemma ll.Sl and Corollary 11.91 are readily generalized to the case with (Bl) 
replaced by (Bl') For £ sufficiently small, the inviscid shock (n^,M^) is weakly stable, with 
all pure imaginary roots simple and lying away from the singularities of A^. In this case we 
obtain a famiy of roots/crossings of the imaginary axis, one for each imaginary root of A*". 
In particular, for gas dynamics, due to rotational invariance (see example ll.2p . we obtain 
families of four crossing eigenvalues A±(efc), with each of A+ and A_ occurring at both 
wave-numbers k and —k. This is not a standard Hopf bifurcation, but a more complicated 
version with 0(2) symmetry, and so we cannot apply directly Theorem 11.111 
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1.6 Discussion and open problems 

We have presented in a simple setting a rigorous mathematical demonstration of a mech- 
anism by which destabilization of hyperbolic surface waves arising in the inviscid shock 
stability problem in the whole space can, at appropriate transverse length scales, lead to 
Hopf bifurcation of a viscous shock in a finite-cross-section duct: specifically, destabilization 
of the effective transverse viscosity (3 investigated in [ZS[ IBSZj , or violation of the refined 
stability condition. This appears to be a fundamentally different mechansim than the hy- 
perbolic ones proposed by Majda et al [MRU IMR21 lAM] via weakly noninear geometric 
optics. 

We point out that as shock parameters cross the inviscid strong instability boundary, gas- 
dynamical shocks undergo one-dimensional instability, or wave-splitting, a more dramatic 
change in front topology than the cellular instabilities we seek to investigate. Thus, cellular 
instability must occur before the strong instability boundary is reached. Experimental 
observations, though not conclusive, indicate that nonetheless it occurs near the strong 
instability boundary [BEj, suggesting that it lies in the region of neutral inviscid stability 
as we have conjectured. 

More, if the transition to cellular instability occurs at low frequencies, it must occur by 
the scenario described, or else the viscous shock would remain stable up to the point of 
wave-splitting. If, on the other hand, it occurs at high frequencies, then as pointed out in 
|BSZ| . then it necessarily involves Hopf bifurcation, by one-dimensional inviscid stability, 
(H2) (satisfied for typical equations of state). Thus, it would appear quite promising to 
search for Hopf bifurcations in the region of neutral inviscid stability, whether of the "low- 
frequency" type studied here or a "high-frequency" type involving unknown mechanisms. 
This would be a very interesting direction for numerical investigations, for example by the 
numerical Evans function techniques of [Brit IBr2[ IBrZ[ IBDG[ IHuZ[ IBHRZj . 

Another interesting direction would be investigation of the stability coefficient /?, both 
numerically and analytically. As pointed out in [BSZj . this is numerically quite well- 
conditioned. One might also consider attempting to carry out an asymptotic analysis near 
the endpoints of the region of neutral stability, at which the imaginary root approaches 
either a branch point of A or else infinity. 

At a technical level, an interesting open problem is to carry out a bifurcation analysis 
in the rotationally symmetric case, for example, for gas dynamics, in which the bifurcation 
associated with crossing X± no longer a standard Hopf bifurcation but a more complicated 
type involving 0(2) symmetry. For a description of Hopf bifurcation with 0(2) symmetry, 
see, for example, [Wj . For a circular cross-section, there is besides axial translation an 
additional continuous group invariance of rotation in the transverse direction, leading to 
the possibility of "spinning" instabilities. These degenerate cases require further analysis at 
the level of the finite-dimensional reduced equations; however, the initial reduction to finite 
dimensions, as carried out in [TZ2| . is essentially the same. Other open problems are to 
extend to detonations, as done for the one-dimensional case in [TZ4j and to treat also real, 
or partial viscosities. As discussed in [TZ2[ [TZ3] . the latter problem involves interesting 
issues involving Lagrangian vs. Eulerian formulations. 



11 



2 Conditional stability analysis 

Nonlinear stability follows quite similarly as in the one-dimensional case [Z4j . decomposing 
behavior into a one-dimensional (averaged in X2) flow driving time-exponentially damped 
transverse modes 

Define the perturbation variable 

(2.1) v{x,t) :=u{x + a{t),t) -u{x) 

for u a solution of (|1.2|) - (|1.3|) . where a is to be specified later. Subtracting the equations 
for u{x + a{t),t) and u{x), we obtain the nonlinear perturbation equation 

2 

(2.2) vt- Lv = ^Nj{v)x^ + dta{ux^ + d^^v), 
where 

2 

(2.3) L := A, - ^ d^^A\x), := dfj{u) 

i=i 

denotes the linearized operator about u and Nj{v) := —{f^{u + v) — f^{u) — df^ {u)v), where, 
so long as \v\fji (hence \v\loc. and |u|ioo) remains bounded, 

(2.4) N^{v) = Oi\v\^), d,N^{v) = 0{\v\\d,v\), d^N^ {v) = 0{\d,\^ + \v\\dlv\). 
2.1 Projector bounds 

Let n„ denote the eigenprojection of L onto its unstable subspace and lies = Id — 
the eigenprojection onto its center stable subspace T,cs- 

Lemma 2.1. Assuming (HO)-(Hl), there is Hj defined in (|2.7p such that 

(2.5) Ujd^ = d^Uj 
for j = u, cs and, for all 1 < p < 00, < r < 4, 

(2-6) - ~ 

Proof. Recalling that L has at most finitely many unstable eigenvalues, we find that Hu 
may be expressed as 

p 



^ Indeed, though we do not do it here, this prescription could be followed quite literally at the nonlinear 
level. 
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where <f)j, j = l,...p are generalized right eigenfunctions of L associated with unstable 
eigenvalues Xj, satisfying the generalized eigenvalue equation (L — Xjy^(f)j = 0, rj > 1, and 
(pj are generalized left eigenfunctions. Noting that L is divergence form, and that Xj ^ 0, 
we may integrate (L — XjY^cjij = over R to obtain A^-^ J (j)jdx = and thus J cpjdx = 0. 
Noting that (pj, and derivatives decay exponentially in xi by separation of variables and 
standard one-dimensional theory \H.e\ IZH[ IMaZlj . we find that (j)j = dx^j with and 
derivatives exponentially decaying in xi, hence 

(2.7) Uuf = Y,'^j{dx4>J)- 

3 

Estimating \diUuf\Lp = \ Y^j di(l)j{(f)jf)\Lp < J2j \di4>j\Lp\(pj\Li\f\Lp < C[/|lp for l/p + 
1/q = 1 and similarly for d^Huf, we obtain the claimed bounds on n„ and n^, from which 
the bounds on lies = Id — II^ and lies = Id — II^ follow immediately. □ 

2.2 Linear estimates 

Let Gcs{x,t;y) := Ilcse^'^6y{x) denote the Green kernel of the linearized solution operator 
on the center stable subspace Scs- Then, we have the following detailed pointwise bounds 
established in [TZ2llMaZl] . 

Proposition 2.2 r [TZ2l[MaZT] l. Assuming (H0)-(H2), (Dl)-D(3), the center stable Green 
function may he decomposed as Gcs = E + G, where 

(2.8) E{x,t;y) = dxMxi)e{yi,t), 

p.) efe.., . E (^-^ (^) - (^)) tto) 

for 2/1 < and symmetrically for yi > 0, G M" constant, and are the eigenvalues of 
df[u±), and 

(2.10) \ j^dlG{-,t-y)f{y)dy\LP <G{l + t-'^)t-"^^hj.\f\L,, 

(2.11) |^9^G',(.,t;y)/(y)dyUp <C(l + rf)t-^(i-^)-^|/|i„ 
for all t > 0, < s < 2, some G > 0, for any 1 < q <p and f £ L'^ n W . 



13 



Proof. As observed in |TZ2| . it is equivalent to establish decomposition 

(2.12) G = Gu + E + G 

for the full Green function G{x,t;y) := e^^5y{x), where 

p 

Gu{x,t;y) := n,e^%(x) = e^* (/>,(x)<^,(y)* 

i=i 

for some constant matrix M G C^^^ denotes the Green kernel of the linearized solution 
operator on S^, and cpj right and left generalized eigenfunctions associated with unstable 
eigenvalues Xj, j = 1, . . . ,p. 

Using separation of variables, moreover, we may decompose G = G'^, where G^ is the 
Green function acting on Fourier modes of wave number k, i.e., G^ = J-~^G6{^ — TTk/M)J^, 
where J- denotes Fourier transform in X2, and ^ Fourier frequency. This reduces the problem 
to that of deriving the asserted bounds separately on the one-dimensional Green function 
G^ and on the complement "^^^oG^, where the difficulty, due to lack of spectral gap, is 
concentrated in the estimation of the one-dimensional part G^. 

The bounds on the one-dimensional Green function G*^ have already been established in 
|Z4| . Proposition 4.2, by essentially the same stationary phase estimates used in |MaZ3] in 
the stable case XI^j = 0; see |TZ21[Z4] for further discussion. The bounds on the complement 
X^fc^^o follow by the straightforward semigroup estimate |e^*/|LP < Ce-''*|/|LP, ry > 0, 
where L denotes the projection of L onto the intersection of its center stable subspace and 
the subspace of functions with transverse Fourier wave numbers ^ 0, which evidently has a 
nonzero spectral gap (7{L) < —2rj < for some > 0; see |TZ2] for related computations. 

□ 

Corollary 2.3 ( [Z4J ). The kernel e satisfies for all t > 

|e,(-,t)|Lp,|et(-,t)|iP <Ct-5(i-Vp), 
|ei,(-,t)|L. <Ct-5(i-Vp)-i/2. 
Proof. Direct computation using definition (j2.9p . □ 



2.3 Reduced equations 

Recalling that dxiU is a stationary solution of the linearized equations ut = Lu, so that 
Ldx^u = 0, or 

G{x,t;y)ux^{yi)dy = e^^Ux^{xi) = dx^u{xi), 



c 



we have, applying Duhamel's principle to ()2.2 
v{x,t) = / G{x,t;y)vo{y)dy 



c 



JC 



Gy{x, t - s; y){N{v) + dv){y, s) dy ds + a{t)dxiu{xi). 
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Defining 

a{t) = - I e{y,t)vo{y)dy 



c 



(2.13) 

+ 11 ey{y,t - s){N{v) + av){y,s)dyds, 



t 

Jc 



following [ZHt IZ4t IMaZ2t [MaZ3j . where e is defined as in (j2.9p . and recalling the decompo- 
sition G = E + Gu + G of (j2.12p , we obtain the reduced equations 



v{x,t) = I {Gu + G){x,t;y)vo{y)dy 

(2.14) 



c 



JC 



{Gu + G)y{x, t - s; y)iN{v) + av){y, s)dy ds, 



and, differentiating (|2.13p with respect to t, and observing that ey{yi,s) — > as s — > 0, as 
the difference of approaching heat kernels, 



a{t) = - I et{y,t)vo{y)dy 
(2.15) K 

+ / eyt{y,t- s){N{v) +av){y,s)dyds. 



JC 



2.4 Nonlinear damping estimate 

Proposition 2.4 ( [MaZ3j ). Assuming (H0)-(H3), let vq € , and suppose that for < 
t <T, the norm of v remains bounded by a sufficiently small constant, for v as in (|2.ip 
and u a solution of (jl.2p - (|1.3p . Then, for some constants 61^2 > 0, for all < t <T, 

(2.16) 11^^(011^2 < Ce''-^'\\v{0)\\l2 +G f e-'-^'-'\\v\l, + \a\^){s) ds. 

Jo 

Proof. Energy estimates identical with those of the one-dimensional proof in [Z4j , using the 
fact that boundary terms in X2 are identically zero due to periodic boundary conditions. □ 

2.5 Proof of nonlinear stability 

Decompose the nonlinear perturbation v as 

(2.17) v{x,t) = w{x,t) + z{x,t), 
where 

(2.18) w := HcsV, z := HuV- 
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Applying lies to ()2.14p and recalling commutator relation ()2.5p . we obtain an equation 

w{x,t) = / G{x,t;y)wo{y)dy 

(2.19) ^ ^ 

- / Gy{x,t - s;y)Ucs{N{v) + dv){y,s)dy ds 
Jo Jc 

for the flow along the center stable manifold, parametrized hy w & Scs. 

Lemma 2.5. Assuming (HO)-(Hl), for v lying initially on the center stable manifold Mcs, 

(2.20) \z\w-^v < C\w\jj2 

for some C > 0, for all 1 < p < oo and < r < 4, so long as \w\ff2 remains sufficiently 
small. 

Proof. By tangency of the center stable manifold to Scs, we have immediately \z\fj2 < 
C\w\'j^2, whence (j2.20p follows by equivalence of norms for finite-dimensional vector spaces, 
applied to the p-dimensional subspace S^. (Alternatively, we may see this by direct com- 
putation using the explicit description of HuV afforded by Lemma l2.1i ) □ 

Proof of Theorem I i.^l Recalling by Theorem 11.11 that solutions remaining for all time in a 
sufficiently small radius neighborhood of the set of translates of u lie in the center stable 
manifold A4cs, we obtain trivially that solutions not originating in Aics must exit Af in 
finite time, verifying the final assertion of orbital instability with respect to perturbations 
not in Aics- 

Consider now a solution v G Aics, or, equivalently, a solution w G T,cs of (|2.19p with 
z = ^cs{w) G Su. Define 

(2.21) C{t):= sup (\w\h2{1 + s)^ + {\w\Loo + \a{s)\){l + s)^). 

0<s<t ^ ' 

We shall establish: 

Claim. For alH > for which a solution exists with C, uniformly bounded by some fixed, 
sufficiently small constant, there holds 

(2.22) C(t) <C2(^0 + C(i)') for := \v^\^r^u^ . 



From this result, provided < 1/4(71, we have that C(t) < 2C2EQ implies (^{t) < 
2C2-E0) and so we may conclude by continuous induction that 

(2.23) C{t) < 2C2E0 

for all t > 0, from which we readily obtain the stated bounds. (By standard short-time 
existence theory, v G exists and C remains continuous so long as C remains bounded by 
some uniform constant, hence (j2.23p is an open condition.) 
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Proof of Claim. By (|2.6p . ItuoliinH^ = incsfoliini/^ < CEq. Likewise, by Lemma [23] 
([221]), ([23]), and Lemma EH for < s < i, 



(2.24) 



|n,,(A^(^;) + av)iy, s)\l2 < CC(t)'(l + s) 



Combining the latter bounds with representations (j2.19p and (j2.15p and applying Propo- 
sition [221 we obtain 



\w{x,t)\LP < I G{x,t;y)wo{y)dy 
ft 



LP 



(2.25) 



+ 



Jc 



Gy{x, t - s; y)Ucs{N(v) + av){y, s)dy ds 



LP 



<EQ{l + ty^2^^'l'^ + CC{tf / {t- sy^^^^{l + s)-"idyds 

Jo 

<C{Eo + C{tf){l + t)-^^^'--p^ 



and, similarly, using Holder's inequality and applying Corollary 12.31 

Ht)\ < jjetiy,t)\\vo{y)\dy 
ft I- 

+ I I \eyt{y,t - s)\\{N{v) +dv){y,s)\dyds 



Jc 



(2.26) 



< |et|L°c|tio|Li + C'C(*)^ / leytlL'^i't - s)\{N{v) + av)\L2{s)ds 

Jo 

<Eoil + t)--2 +Catf I (t-s)'i{l + s)-Us 

Jo 

<c{Eo+at)^)ii+t)~"^- 



By Lemma 
(2.27) 



In particular, \z\i2{t) < CC(t)^(l + t) 2. Applying Proposition 12.41 and using (j2.25p and 
(j2.26p . we thus obtain 



(2.28) 



\w\H2{t)<C{Eo + atf)il + t)—^. 



Combining (j2.25p . (|2.26p . and (|2.28p . we obtain (|2.22p as claimed. As discussed earlier, 
from (j2.22p . we obtain by continuous induction (]2.23p . or C < 2C2|fo|Lin_ff2, whereupon the 
claimed bounds on \v\lp and \v\jj2 follow by (|2.25p and ()2.28p . and on [d;| by ()2.26p . Finally, 
a computation parallel to (IZlHIl (see, e.g., |MaZ3[[Z2] ) yields \a{t)\ < C{Eq + C,{tf ), from 
which we obtain the last remaining bound on |a(t)|. □ 



17 



References 



[AT] G. Abouseif and T.Y. Toong, Theory of unstable one- dimensional detonations, Combust. Flame 

45 (1982) 67-94. 

[AGJ] J. Alexander, R. Gardner and C.K.R.T. Jones, A topological invariant arising in the analysis 

of traveling waves, J. Reinc Angew. Math. 410 (1990) 167-212. 

[AIT] R. L. Alpert and T.Y. Toong, Periodicity in exothermic hypersonic flows about blunt projectiles, 

Acta Astron. 17 (1972) 538-560. ' 

[AM] M. Artola and A. Majda, Nonlinear kink modes for supersonic vortex sheets, Phys. Fluids A 1 

(1989), no. 3, 583-596. 

[BRSZ] S. Benzoni-Gavage, F. Rousset, D. Serre, and K. Zumbrun, Generic types and transitions 
in hyperbolic initial-boundary-value problems, Proc. Roy. Soc. Edinburgh Sect. A 132 (2002) 
1073-1104. 

[BE] A.A. Barmin and S.A. Egorushkin, Stability of shock waves, Adv. Mech. 15 (1992) No. 1-2, 

3-37. 

[BHRZ] B. Barker, J. Humpherys, , K. Rudd, and K. Zumbrun, Stability of viscous shocks in isentropic 
gas dynamics, Comm. Math. Phys. 281 (2008), no. 1, 231-249. 

[BSZ] S. Benzoni-Gavage, D. Serre, and K. Zumbrun, Transition to instability of planar viscous shock 

fronts, to appear, Z.A.A. 

[BMR] A. Bourlioux, A. Majda, and V. Roytburd, Theoretical and numerical structure for unstable 

one-dimensional detonations. SIAM J. Appl. Math. 51 (1991) 303-343. 

[BDG] T.J. Bridges, G. Berks, and G. Gottwald, Stability and instability of solitary waves of the 

fifth-order KdV equation: a numerical framework. Phys. D 172 (2002), no. 1-4, 190-216. 

[Brl] L. Brin, Numerical testing of the stability of viscous shock waves. Doctoral thesis, Indiana 

University (1998). 

[Br2] L. Q. Brin, Numerical testing of the stability of viscous shock waves. Math. Comp. 70 (2001) 

235, 1071-1088. 

[BrZ] L. Brin and K. Zumbrun, Analytically varying eigenvectors and the stability of viscous shock 

waves. Seventh Workshop on Partial Differential Equations, Part I (Rio de Janeiro, 2001). Mat. 
Contemp. 22 (2002), 19-32. 

[Erl] J. J. Erpenbeck, Stability of step shocks. Phys. Fluids 5 (1962) no. 10, 1181-1187. 

[Fl] W. Fickett, Stability of the square wave detonation in a model system. Physica 16D (1985) 

358-370. 

[F2] W. Fickett, Detonation in miniature, 133-182, in The mathematics of combustion. Frontiers in 

App. Math. (1985) SIAM, Philadelphia ISBN: 0-89871-053-7. 

[FW] Fickett and Wood, Flow calculations for pulsating one- dimensional detonations. Phys. Fluids 

9 (1966) 903-916. 



18 



[Fo] G.R. Fowles, On the evolutionary condition for stationary plane waves in inert and reactive 

substances, in Shock induced transitions and phase structures in general media, 93-110, IMA 
Vol. Math. Appl., 52, Springer, New York, 1993. 

[FT] H. Freistiihler and Y. Trakhinin, On viscous and inviscid stability of magnetohydrodynamic 

shock waves, to appear, Physica D. 

[GZ] R. Gardner and K. Zumbrun, The Gap Lemma and geometric criteria for instability of viscous 

shock profiles. Comm. Pure Appl. Math. 51 (1998), no. 7, 797-855. 

[He] D. Henry, Geom,etric theory of semilinear parabolic equations. Lecture Notes in Mathematics, 

Springer Vcrlag, Berlin (1981), iv + 348 pp. 

[HuZ] J. Humpherys and K. Zumbrun, An efficient shooting algorithm for Evans function calculations 

in large systems, Phys. D 220 (2006), no. 2, 116-126. 

[KS] A.R. Kasimov and D.S. Stewart, Spinning instability of gaseous detonations. J. Fluid Mech. 

466 (2002), 179-203. 

[Ml] A. Majda, The stability of multi- dimensional shock fronts - a new problem for linear hyperbolic 

equations. Mem. Amer. Math. See. 275 (1983). 

[M2] A. Majda, The existence of multi- dimensional shock fronts. Mem. Amer. Math. Soc. 281 (1983). 

[M3] A. Majda, Compressible fluid flow and systems of conservation laws in several space variables. 

Springer- Verlag, New York (1984), viii-|- 159 pp. 

[M4] A. Majda, Compressible fluid flow and systems of conservation laws in several space variables. 

Springer- Verlag, New York (1984), viii+ 159 pp. 

[MRl] A. Majda and R. Rosales, A theory for the spontaneous formation of Mach stems in reacting 

shock fronts, I. The basic perturbation analysis, SIAM J. Appl. Math. 43 (1983), no. 6, 1310- 
1334. 

[MR2] A. Majda and R. Rosales, A theory for spontaneous Mach-stem formation in reacting shock 

fronts. II. Steady-wave bifurcations and the evidence for breakdown. Stud. Appl. Math. 71 
(1984), no. 2, 117-148. 

[MaZl] C. Mascia and K. Zumbrun, Pointwise Green's function bounds and stability of relaxation 

shocks. Indiana Univ. Math. J. 51 (2002), no. 4, 773-904. 

[MaZ2] C. Mascia and K. Zumbrun, Stability of small-amplitude shock profiles of symmetric hyperbolic- 
parabolic systems. Comm. Pure Appl. Math. 57 (2004), no. 7, 841-876. 

[MaZ3] C. Mascia and K. Zumbrun, Pointwise Green function bounds for shock profiles of systems with 
real viscosity. Arch. Ration. Mech. Anal. 169 (2003), no. 3, 177-263. 

[MT] U.B. McVey and T.Y. Toong, Mechanism of instabilities in exothermic blunt-body flows, Com- 

bus. Sci. Tech. 3 (1971) 63-76. 

[Me] G. Metivier, Stability of multidimensional shocks. Advances in the theory of shock waves, 25- 

103, Progr. Nonlinear Differential Equations Appl., 47, Birkhauser Boston, Boston, MA, 2001. 



[MeZ] G.Metivier-K. Zumbrun, Symmetrizers and continuity of stable subspaces for 

hyperbolic boundary value problems. Discrete Contin. Dyn. Syst. 11 (2004), no. 1, 205-220 



19 



[R] M. Raoofi, asymptotic behavior of perturbed viscous shock profiles, J. Hyperbolic Differ. 

Equ. 2 (2005), no. 3, 595-644. 

[Sel] D. Serre, La transition vers I'instabilite pour les ondes de chocs multi-dimensionnelles, Trans. 

Amer. Math. Soc. 353 (2001) 5071-5093. 

[Se2] D. Serre, Systems of conservation laws. 1. Hyperbolicity, entropies, shock waves, Translated 

from the 1996 French original by I. N. Sneddon. Cambridge University Press, Cambridge, 
1999. xxii+263 pp. ISBN: 0-521-58233-4. 

[Se3] D. Serre, Systems of conservation laws. 2. Geometric structures, oscillations, and initial- 

boundary value problems. Translated from the 1996 French original by I. N. Sneddon. Cambridge 
University Press, Cambridge, 2000. xii-|-269 pp. ISBN: 0-521-63330-3. 

[TZ2] B. Texier and K. Zumbrun, Galloping instability of viscous shock waves, Physica D. 237 (2008) 

1553-1601. 

[TZ3] B. Texier and K. Zumbrun, Hopf bifurcation of viscous shock waves in gas dynamics and MHD, 

Arch. Ration. Mech. Anal. 190 (2008) 107-140. 

[TZ4] B. Texier and K. Zumbrun, Transition to longitudinal instability of detonation waves is gener- 

ically associated with Hopf bifurcation to time-periodic galloping solutions, preprint (2008). 

[T] Y. Trakhinin, A complete 2D stability analysis of fast MHD shocks in an ideal gas, Comm. 

Math. Phys. 236 (2003) 65-92. 

[W] C. Wulff, Theory of Meandering and Drifting Spiral Waves in Reaction-Diffusion Systems, 

Dissertation, Frcic Universitat Berlin (1996). 

[Zl] K. Zumbrun, Multidimensional stability of planar viscous shock waves. Advances in the theory 

of shock waves, 307-516, Progr. Nonlinear Differential Equations AppL, 47, Birkhauser Boston, 
Boston, MA, 2001. 

[Z2] K. Zumbrun, Stability of large- amplitude shock waves of compressible Navier-Stokes equations, 

with an appendix by Helge Kristian Jenssen and Gregory Lyng, in Handbook of mathematical 
fluid dynamics. Vol. HI, 311-533, North-Holland, Amsterdam, (2004). 

[Z3] K. Zumbrun, Planar stability criteria for viscous shock waves of systems with real viscosity, 

in Hyperbolic Systems of Balance Laws, CIME School lectures notes, P. Marcati ed.. Lecture 
Note in Mathematics 1911, Springer (2004). 

[Z4] K. Zumbrun, Gonditional stability of unstable viscous shocks, preprint (2008). 

[ZH] K. Zumbrun and P. Howard, Pointwise semigroup methods and stability of viscous shock waves. 

Indiana Mathematics Journal V47 (1998), 741-871; Errata, Indiana Univ. Math. J. 51 (2002), 
no. 4, 1017-1021. 

[ZS] K. Zumbrun and D. Serre, Viscous and inviscid stability of multidimensional planar shock 

fronts, Indiana Univ. Math. J. 48 (1999) 937-992. 



20 



